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Abstract 

We derive dispersion relations for K tttt decay, using the Lehmann- 
Symanzik-Zimmermann formalism, which allows the analytic continuation 
of the amplitudes with respect to the momenta of the external particles. No 
off-shell extrapolation of the field operators is assumed. We obtain general- 
ized Omnes representations, which incorporate the vrvr and irK S'-wave phase 
shifts in the elastic region of the direct and crossed channels, according to 
Watson theorem. The contribution of the inelastic final-state and initial- 
state interactions is parametrized by the technique of conformal mappings. 
We compare our results with previous dispersive treatments and indicate 
how the formalism can be combined with lattice calculations to yield physi- 
cal predictions. 
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1 Introduction 



The weak decay K txt: has been a continuous challenge for the theo- 
retical investigations. Chiral perturbation theory (ChPT) was extensively 
used [0, but the large number of unknown counterterms and renormaliza- 
tion constants render the numerical predictions difficult beyond the leading 
order. Most lattice calculations (see and references therein) simulate ma- 
trix elements of the type (7r|0j|i^'), related to the physical matrix elements 
(7r7r|Oi|A') by lowest order ChPT ||^. In this procedure the higher order final 
state interactions are completely missing, while it is expected that they play 
an important role for the A/ = 1/2 rule and the CP- violating ratio e'/e. 
The finite-volume techniques developed in Q can take into account FSI, but 
they are numerically very demanding [Q. In a combined approach proposed 
recently, the results obtained by lattice simulations at unphysical points are 
extrapolated to the physical configuration by using calculations to NLO in 
ChPT 

An alternative way to connect the on-shell amplitude to lattice results 
at unphysical points and to spectral functions measured experimentally is 
based on dispersion relations. This formalism was used some time ago for 
the CP-conserving amplitudes in order to explain the A/ = 1/2 rule 0, 
and more recently in p], Q for evaluating the effects of final state interac- 
tions upon e'/e. The last works use an Omnes representation |10| for the 
decay amplitude, written by analogy with the case of the scalar form fac- 



tor. This approach was investigated further in Refs. ITll-llT^l, where some 
critical remarks about the method were advanced. An alternative dispersive 
framework for K tttt decay was proposed in [|14|, by assuming that the 
weak hamiltonian carries a non-zero momentum. Then the matrix element 
of the decay becomes equivalent with the nK elastic scattering amplitude, 
for which Mandelstam representation is assumed. 

In the references mentioned above, the dispersion relations for the weak 
decay were written by using the analogy with the familiar cases of the form 
factors or the scattering amplitudes. However, in the weak decay a contin- 
uation in the external momenta is necessary in order to obtain a dispersion 
relation. As a proof of the dispersion relations in this case is missing, their 
meaning was not always clear and led to some confusion. The dispersive 
variable was interpreted either as the mass or the momentum of an off-shell 
particle. The clarification of this point is possible only by a systematic deriva- 
tion in the frame of a field theoretic formalism. In the present work we ad- 
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dress this problem, by performing the continuation in the external momenta 
with the Lehmann-Symanzik-Zimmermann (LSZ) formalism |T^. Our main 
result is a general Omnes representation for the K —>■ tttc amplitudes, includ- 
ing final and initial state interactions in both the direct {K — >■ tttt) and the 
crossed (tt —>■ ttK) channels. The derivation clarifies the significance of the 
dispersive variables, allowing to make contact with lattice calculations done 
at unphysical points. 

In the next section we present the derivation of the dispersion relations, 
using LSZ reduction and hadronic unitarity. We follow to some extent the 
dispersive treatment of B ^ tttc decay considered in [O, Il8|. However, the 



different masses of the decaying particles in the two processes require specific 
treatments. In section 3, we derive a generalized Omnes representation for 
K ^ TTTT decay by solving the inhomogeneous Hilbert problem [|19], ^ 



m 



the direct and the crossed channels. In Section 4, we compare our results 
with the dispersion relations considered previously in the literature, and show 
how to combine them with lattice calculations in order to predict the physical 
amplitude. 



2 LSZ reduction and dispersion relations 

We consider decay amplitudes Aj of definite isospin, 1 = 0,2, defined as 

Aj = {{Tr{h) Tr{k2))i ; out \nUO)\K{p) ; in) , (1) 

where the "in" and "out " states are defined with respect to the strong in- 
teractions and Hw is the weak effective hamiltonian density pi 



HJx) 



k=u,c 



(2) 



Ci(/i)Of(x,/i) + C2(/x)02^(a;,/i) 



J2 Cj{fi)Oj{x,fi) 

j=3,...,8 



Here Oj are local AS = 1, AB = operators and Cj the corresponding Wil- 
son coefficients, which take into account perturbatively the strong dynamics 
at distances shorter than We assume that a factor i was included in the 
definition of the operators, so that the amplitudes (|T]) satisfy time reversal 
invariance, up to the complex coefficients in (^). 
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For our purpose it is more convenient to start from the 5'-matrix element 



Si = (vr(A;i) 7r(A;2); out \K{p); in) , (3) 

where the transition from the "in" to the "out" states is achieved by both the 
strong and weak interactions. By expanding the S*— matrix to first order in 
the weak interactions one obtains the expression (0) of the decay amphtude. 
Alternatively, by applying the LSZ reduction [jl6[ to the K meson in Eq. 
the decay amplitude (|l|) is expressed as 

Aj = -l=(7r(A;i) 7r(A;2); out | riK{0) |0) , (4) 

where riK{x) = ICx4>k{x) denotes the source operator {JCx is the Klein-Gordon 
operator and (px the interpolating field of the kaon). In a Lagrangian theory 
the source operator has the formal expression 

d(pK od^cpK 

i.e. it has contributions from both the strong and weak parts of the inter- 
action Lagrangian. In what follows we do not need the explicit expressions 
of the sources, but only the significance of the matrix elements involving 
them. We stress also that throughout the derivation the sources are on-shell 
operators, defined in terms of the physical interpolation fields. 

The matrix element depends on the momenta ki and k2 of the two 
pions. We shall consider it as a function of the invariant variables s = 
(fci + ^2)^, t = kf, u = k\. The physical amplitude corresponds to the values 
s = mj^, t = and u = m^. The extrapolation to arbitrary external 
momenta can be achieved by the LSZ reduction formalism ||T6|. We remark 
that Eq. @) is similar to the definition of the electromagnetic form factor 
of the pion, where tik is replaced by the electromagnetic current J^. We 
can apply therefore the standard methods used in deriving the dispersion 



relations for the pion form factor Making the LSZ reduction of one final 
pion in Eq. (H), we obtain 

Ai{s,t) = -^^= J da;e*^^(xo) (7r(A;2)| 77^(0)] |0) , (6) 

where riT,{x) is the source of the reduced pion. We left aside the so-called 
"degenerate terms" which are polynomial of the Lorentz invariant variables 
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Then Eq. (^) defines a function liolomorpliic for those values of the 
external squared momenta s and t for which the integral is convergent (for 
the unreduced pion we take the physical value u = ml). Due to the presence 
of 9{xq), the integral upon xq converges in the upper half of the fcio complex 
plane, Imfcio > 0. The causality property of the commutator restricts the 
integral upon the spatial variables to |x| < |xo| [Q. We choose the particular 
Lorentz frame where the unreduced pion is at rest (k2 = 0), when kiQ = 
{s-t-ml)/2m^ and = [s - {Vi + m^f][s - {Vi - m^)^]/{2m^y. Then 
the integral in (|]) represents a function of s and t, analytic for complex 
values of these variables, with possible discontinuities along the real axes. 
The rigorous proof of the analyticity in the external masses is actually a 
difficult problem and requires a more detailed analysis |2^. Here we do not 



attempt to give a proof, but only use the LSZ representation to understand 
the meaning of the dispersive variables and to read off the contributions to 
the spectral functions appearing in the dispersion relation. 

The discontinuity across the real axis is obtained formally from the expres- 
sion (ID by replacing i9{xo) by 1/2, inserting a complete set of intermediate 
states in the commutator [riT^{x),'r]K{0)] and using translational invariance 



22| . The two terms in the commutator allow us to decompose the spectral 

a = a^ + at, (7) 



function as 
where 



^ J26{h + k2-Pn){Ak2)\vA0)\n){n\r]Ki0m , (8) 



and 



= ^-^=_E^(^i +Pn)(vr(A;2)|r7i,(0)|n)(n|r^.(0)|0) . (9) 

In these relations the summation is over intermediate states consisting of 
physical particles, with an implicit integration upon their momenta. By 
the subscripts s (t), we anticipate the fact that ag receives contributions 
from the s— channel {K — »• tttt) and at from the t— channel (tt —>■ Ktt). 
Accordingly, we can write the amplitude as a sum of two terms, Al and Aj., 
obtained by a dispersion representation involving the spectral function ag 
and at, respectively. In order to evaluate the spectral functions, we recall 
that the sources contain contributions from both the strong and the weak 
interactions, the last ones being treated to first order. 
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Let us consider first the spectral function as defined in (|^). As discussed in 
Ref. Il^], the intermediate states n which contribute to the unitarity sum are 



generated by either the weak or the strong part of the source rjK, undergoing 
a rescattering to the final tttt state by a strong (weak) process, respectively. 
The first contributions represent the so-called "final state interactions" (FSI), 
while the second are usually interpreted as "initial state interactions" (ISI). 
The lowest intermediate state contributing to FSI consists of two-pions, which 
produces the branch-point s = ^m^, while for the ISI the lowest intermediate 
state is the pair K^-k, responsible for the threshold s = {rriK* + m.,,)'^. In 
order to write the specific contributions, we recall that in the LSZ formalism 
the matrix elements of the sources represent, up to kinematical factors, the 
physical decay or scattering amplitudes |]22|. Thus, according to Eq. (^, 



(7r7r|?7;^(0)|0) is the amplitude of the weak decay of K into two pions (here 
only the weak part of the source contributes), while (7r(A;2)|^?7r(0)|n) is the 
amplitude of either the strong or the weak n — > tttt transition, depending on 
which part of the source is considered. 

A remarkable property of cr^ is that it depends only on s, being indepen- 
dent of the variable t fl^, [^. This can be easily seen by recalling that the 



intermediate states \n) consist of physical particles. By choosing the cm. 
system, where = s is the total energy squared, we see that the matrix el- 
ements in (|^) depend only on s and the physical masses. In the two-particle 
approximation of the unitarity sum, the integral in (§) can be performed 
exactly |T^. According to the discussion above we can write 

as{s) = cTpsiis) + (Tisi{s) , (10) 

where the first contributions to each term are 

(TFSiis) = e{s - AttiDM^^^^^Ak-..^ + e{s - Aml)M*ij,^^^Aj,^jij, + . . . , 
cTisiis) = e{s - {niK* + m^y)A*^,^_,^^MK-.K*As) + ■■■ ■ (H) 

In the above relations M^^^^^ and Mj^j^_^^^ denote on-shell S-wave strong 
scattering amplitudes at cm. energy squared s, and Ak^-kw {Mk^k'-k), etc., 
are weak (strong) decay amplitudes. The amplitude Al can be recovered 
from the discontinuity by means of a dispersion integral. Neglecting for the 
moment possible subtractions and polynomials in the Mandelstam variables, 
we have 
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We consider now the spectral function at defined in (|^). The weak part 
of the source rj^^ is responsible for the FSI in the t-channel, with the lowest 
branch-point at t = {rriK + '^tt)^- The strong part of the source r^jr generates 
the ISI in the t-channel, with the lowest branch-point t = 9m^. As above, it 
is easy to show that at does not depend on s and can be written as 

Mt) = ^Fsiit) + aisiit) , (13) 

where, according to the above discussion, the lowest terms are 

^Fsi{t) = 9{t - {niK + m^)^)A^*^_^^^A^_^i<- + . . . , 

&isi{t) = e{t-'dml)N:^^^A^^^^K + ... . (14) 

Here N^^k^-kK is the iiK S-wave scattering amplitude at cm. energy squared 
equal to t, and A^^^^^ (A^^^3^) etc., are weak (strong) decay amplitudes. 
Neglecting again possible subtractions, we write the amplitude in terms 
of its discontinuity in the t-channel by a dispersion integral 

The total amplitude Ai{s,t) is then expressed as the sum 

Aj{s,t) = Aiis) + Alit) , (16) 

the physical amplitude being obtained for s = m\ and t = m^. We note in 
particular that Af(m^) is a real number, since the point t = is situated 
below the cuts in the dispersion relation ([T5|). 

The significance of the variables s and t is clear from the above discussion: 
s is defined in terms of the pion momenta as s = {ki+k2)^ and t is equal to the 
external momentum squared of one pion, t = k\. Therefore, it represents the 
mass squared of one external pion. We recall that in the above formalism no 
off-shell extrapolation was assumed, the sources entering the matrix elements 
being on-shell operators. 



3 Omnes representations 

It is convenient to write the above dispersion relations in terms of the phases 
of the rescattering amplitudes in the elastic region, according to Watson theo- 



rem [24|. To illustrate the method, we consider first the amplitude ^/(s, m^) 
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as a function of s at physical t = m^. The general case of Aj{s,t) will be 
treated in subsection pl2| . This generalization is useful in order to incorporate 
information available on the decay amplitude at nonphysical pion masses. 



3.1 Amplitude ^4/(5, m^) 



In this mentioned above, the last term in Eq. ([Iq) is a real constant. 

Denoting A± = Aj{s±ie, ml), we write the unitarity relation in the s-channel 
as 



A^-A^ 



2i 



where 



Mi(s) 



1 



2i 



(17) 



(18) 



is the S-wave nn scattering amplitude of isospin /. In Eq. (|T^) crj„ denotes 
the sum of the inelastic FSI and the ISI spectral functions (we take Sj„ equal 
to the FSI inelastic branch-point 4m|^, which is lower than the ISI branch- 
point {rriK* + ^7r)^)- In the r.h.s. of Eq. ( pTj) we note the presence of the 
amplitude A^ = Af{s + ie,m1), due to the fact that the intermediate two 
pions in the unitarity sum are physical particles, as we mentioned above. 
The relation (|17D can be written as an inhomogeneous Hilbert equation 

PI 



A+{1 - 2iM*) -A^ = 2i e{s - Sin) cT^nis) , s > 4m^ 



(19) 



We shall construct the solution by imposing time reversal invariance, which 
implies that the amplitudes satisfy the reality condition Af{s*) = A}{s) and 
the discontinuity across the cut is equal to the imaginary part [^. Using 
the expression (|T8|) we obtain from (|19|), for s > 4m^, 



ImAi cos 5q — Re A J sin 5q = 6{s — Si 



5^1 



^ 2Re[a,„e" ^o 
1 + 



(20) 



We define now the Omnes function 



^lis) 



exp 



S - So 



TT 



5^(s') ds' 



Ami 



{s' - So){s' - s) 



(21) 



assuming that one subtraction is sufficient. The boundary values of f2/(s) 
satisfy the relations fi/(s±ie) = exp(±z5g)|n/(s)| (we recall that the modulus 
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|f2/(s)| is obtained from PT| ) by taking the principal value of the integral). 
Then Eq. (EH) can be written as 



Im 



^■'-"'^TT^^kMT- '''' 



a,(s) 

We define now the function Gi{s) through the relation 

Ms,ml) = nj{s)Gj{s), (23) 

and express it by a dispersion relation in terms of its imaginary part given 
in dl): 

Here Pi{s) is a polynomial and, for convenience, we wrote the integral with 
one subtraction. The subtractions are actually not relevant in our method, 
since we shall parametrize the function Gi{s) in a different way, using the 
technique of conformal mappings^. Namely, since by construction Gj{s) is 
analytic in the s— plane cut for s > s^, we consider the variable 



4s) = I ^ ■—= , (25) 



which maps the s-plane cut along the real axis for s > Sin onto the disk z < 1 
of the plane z = z{s). Actually, the mapping of the s-plane onto the unit 
disk is not unique . For further convenience we choosed the mapping such 
that zim^) = 0. Now we expand Gi{s) in powers the variable z 

Gj{s) = j:ai!n<s)r, (26) 

n 

where al^^ are real numbers. This series converges in the whole disk \z\ < 1, 
i.e, in the whole s-plane cut along s > Sj„, in particular at the physical point 
s = m|-. 



^The method of conformal mappings was proposed in particle physics a long time 
ago psf . In a context similar to the present one, the method was applied to the pion 
electromagnetic form factor in Ref. [p6[. 
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Inserting into (^) we obtain a representation of the amplitude 

Ms,ml)=nj{s)Y.^i!\<s)r, (27) 

n 

in terms of the known S*— wave tttt phase shifts entering the Omnes function 
and the real Taylor coefficients a'^\ In Section 5 we shall discuss how 
these coefficients can be determined by using lattice results at unphysical 
values of s. The physical amplitude is obtained from (^) by setting s = m\. 



3.2 Amplitude ^/(s,t) 

We consider now the amplitude Ai{s, t) for arbitrary arguments. It turns out 
that the elastic unitarity for t) can not be immediately solved by means 
of an Omnes representation, as in the above treatment of y4/(s,m^). Indeed, 
from the relations (|T^) and the discontinuity of Ai{s,t) across the 

cut along s > 4m^ at fixed t is 

Aj(s + ie,t) — Aj(s — ie,t) , . sn /,/ n / n 

-^^—^ — = M;{s)Aj{s + te,ml) + e{s - Sin)cTin{s) . (28) 

We note the presence, in the r.h.s., of the amplitude Ai{s + ie^m^), due to 
the fact that the intermediate states in the unitarity sum consist of physical 
particles. Therefore the functions which appear on the two sides of (^) are 
different, and this relation can not be written as a Hilbert boundary value 
equation. 

This difficulty can be circumvented if we treat separately the functions 
Al[s) and Al{t) defined in Eqs. ( p!2|) and (|15D, respectively. Denoting As^± = 
Al{s ± ie) we obtain from Eqs. ([Ty-|T^ 

= 9{s - 4ml)M^{s)Aj{s + ze, ml) + 9{s - s,n) cr.„(s) . (29) 

2i 

In order to bring this relation to a form convenient for the Muskhelishvili- 
Omnes technique, we express, according to Eq. (|T6|), 



Aj{s + le, ml) = + Al{ml) . (30) 
Then Eq. (^) becomes: 

^,+ (1 - 2iM*) - = 2i [e{s - Aml)M*A[{ml) + e{s - s,„) a,„(s)] . (31) 
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This equation is similar to Eq. (|T^, excepted for an additional term in the 
r.h.s., which contributes to the imaginary part above the elastic threshold 
4m^. Therefore, As{s) will be of the form (^), with the function Gi{s) given 
by a dispersion relation similar to (|2^, containing in addition the term 



S- So 



ds' 



Re[M*{s')e<^'">] 



Ami 



i + 4\ni{s')\{s' -so){s' -s) 



(32) 



where we took into account that Alim, 



is a real constant, as mentioned 
below Eq. (|16|). 

It is convenient to separate in this integral the contribution of the inelastic 
region s > Sm, combining it with the contribution of the inelastic term atn 
and expanding them in powers of the variable z defined in (p5|). Therefore, 
we express Al{s) as 



Aii 



njis) 



with 



Ms) 



s- So 



vr 



sin Sq{s') ds' 
\ni{s')\{s' - so){s' - s) 



(33) 



(34) 



sinSn. 



where we took along the elastic region ?7o(s) = 1 and Mj 

It is convenient to choose the subtraction point Sq = m'^ in both the 
expression ( PT] ) of the Omnes function and the definition ( P^ of the function 
fi{s). This implies f2/(m^) = 1 and //(m^) = 0. By recalling also that the 
conformal mapping ( P^D was defined such as z{m'^) = 0, it follows that the 
amplitude Al{s) given by Eq. (|33|) is normalized as 



Aiiml) = cH^ . 



(35) 



We consider now the second term in the decomposition ([16|) , namely the 
amplitude Al{t). In order to obtain an Omnes representation, we must work 
with amplitudes Aj{s,t) of definite isospin in the t channel, vr — > irK. By 
crossing symmetry we can write 



Ai{s,t) 



E 

2'2 



CijAj{s,t) 



0,2 



(36) 
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where the matrix Cjj is known from the elastic ttK scattering [^. Each 
amphtude Aj{s,t) admits a decomposition similar to ([T6| ) 



(37) 



We consider the amplitude A^{t) and define At^± = ^/(t ± i^)- Then the 



unitarity relation in the t-channel, given by Eqs. (|l3| - p!5| ), can be written as 



2i 



e{t-{mK + m^y)N*j{t)Aj{mi^,t + ie) + e{t-t 

in ) ^in 

(t), (38) 



where 



2i 



(39) 



denotes the S'-wave nK scattering amphtudes at cm. energy squared equal 
to t, and crin{t) the contribution of the inelastic FSI and ISI t-channels. We 
take tin equal to the ISI branch-point 9m^, which is lower than the inelastic 
FSI branch-point {niK + Tnr,Y. 

Since the intermediate tiK state in the unitarity sum (P) consists of phys- 
ical particles, in the r.h.s. of ( pHf ) contributes the amplitude Aj{m\,t + ie). 
By expressing Aj{m\,t + ie) according to ( ^Tf ) we obtain 

A,+(i-2ziV})- A,_ 



= 2i [e{t - [m^ + mKfWjAi{m\) 
+ Q{t-Un)^in{t)\. (40) 



The solution of this equation can be obtained following the procedure applied 
to the function A\{s). We introduce the Omnes function 



rij(t) = exp 



TT 



5i((t'^ 



(m,r+mK)^ 



{t'-u){t'-t) 



dt' 



(41) 



and express the ratio yl/(t)/f2j(t) through a dispersion relation in terms of 
its imaginary part calculated from (pOl). Then Al if) can be written as 



+ 



t-to 



dt'- 



t-u 



vr 



dt'- 



Re[Ai{m\)N}{t')e^'^'o' 



i + fii \nj{t')\{t' -to){t' -t) 



Re [ainity'o 



l + f]i\nj{t')\{t'-to){t'-t) 



(42) 
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We recall that the quantities A'^{mj^) entering this relation are complex 
numbers. 

We further separate in Eq. (^) the contribution of the elastic part of the 
cut, and take into account the higher singularities by means of a conformal 
mapping. Namely, we define the variable 



w{t) 



tin v^in ^ 



tin ytin t 



(43) 



which maps the t-plane cut for t > tin onto the disk \w\ < 1, such that 
w(m^) = 0, and expand the polynomial and the inelastic part of Eq. (^) in 
powers of this variable: 

t-to7^t> 2 Re{[Ai{mj,)N}it') + aUt')]e''o} 



Pj + 



i + vi \nj{t')\it'-to){t'-t) 

''in 

the coefficients c^*^-* being real. Then Eq. ( ^21) becomes 

n 

where we defined 

t-to 



(44) 
(45) 



9j{t) 



sin 5q dt' 



\nj{t')\{t'-to){f-t) 



(46) 



We took into account the fact that in the elastic region fj^ = 1 and NjC = 
sin Sq. 

It is convenient to choose the subtraction point to = in both the 
Omnes function (^) and the definition (^61) of gj{t). This means that 
Qj{m^) = 1 and gj{m'^) = 0. Recalling also that the conformal variable 
w(t) was defined in ( ^31) such that w(m^) = 0, we obtain from (^5]) 

(47) 



~SJ) 
Co 



Collecting Eqs. (ITg,|3^,|36|) and (|45D , we express the amplitude Ai{s,t) 



as 



Ai{s,t) 



J 



n 

ReAi{ml)gj{t) + Y.~c'^;!\^it)\ 

n 
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(48) 



where the functions fi{s) and gj{t) are defined in Eqs. ( Pl] ) and re- 
spectively (we recall that Sj„ = 4m|. and tj„ = 9m^). In Eq. (^8]) we must 
insert, according to (p!6| , |36|j37| ) and (0), 



(49) 



Also, using the crossing relation (|36D, we express the quantity ReA^(m|-) 
entering (|^) as 



Reif(m^) 



E ^7LReAf(m 



K) ? 



J 



L=0,2 



1 3 
2 '2 



(50) 



where ReAf'(m|^) is obtained from (pB]). The relations (^8D-(^) provide a 
system of coupled equations, which express each amplitude Af{s, t) (/ = 0, 2) 
in terms of vrvr and vri^" S'-wave phase shifts and the real coefficients c^^-* and 
c(^)(J= 1/2, 3/2). 

At fixed t = the amplitude (H3) takes the simple form 



Ai(s,ml) = nj(s) 



T.CiA'^ {fi{s) + ni\s)} + Y.^^Ms)T 



I J 



(51) 



where we introduced the last term of Eq. 



inside the brackets, and expressed A[( 



m 



), i.e. A^it) evaluated for 
according to Eq. (|49|) . 



t = 

It is easy to verify that the function multiplying the Omnes factor in the 
relation (|51|) is real for s < Sj„. Indeed, the coefficients c^^^ and c^f^ are real, 
and for values of s below the inelastic threshold the variable z{s) is real. The 
only terms having an imaginary part for s < Sj„ are the function fi{s) and 
the Omnes function Qi{s). But it is easy to check, using (^) and (0), that 
their imaginary parts compensate in Eq. (pTf): 



Im [fj{s)+nj\s)] 



sin (5q(s) 

\ni(s)\ 



sin (5q(s) 



0. 



(52) 



Therefore, the first term in the r.h.s. of (^) is real along the elastic region 
and can be included in the expansion in powers of the variable z{s). This 
shows that the general representation (^8]) reduces, when t = m^, to the 
Omnes representation (^) derived in the previous subsection. For arbitrary 
values of t, however, the amplitude ( ^HD has additional cuts in the elastic 
region of both the s and t channels. 
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The physical amphtude is obtained from (^) for s = m\- and t = rn^. 
With our normahzation, it depends on the S'-wave phase shift of vrvr scattering 
and the Taylor coefficients c^^^ and c^^^ (the S'-wave phase shift of the tiK 
scattering contribute only indirectly, through these coefficients). As proved 
in by the conformal mapping the rate of convergence is improved, so 
we expect to obtain an accurate representation at low energies with a small 
number of terms in the expansions. In the next section we shall discuss how 
to determine the coefficients c^^-* and c^^^ using lattice results at unphysical 
points. 

We end this section with two remarks: first we notice that in the above 
derivation the symmetry between the two final pions is not explicit. This 
symmetry can be easily imposed by writing down a dispersion relation sym- 
metrical with respect to the interchange of t and u. Namely, instead of Eq. 
( [TB| ) we have, more generally 

A,{s, t, u) = Alis) + \ [Alit) + Ai{u)\ , (53) 

where Al^{u) satisfies a dispersion relation similar to (^). 

The second remark concerns the starting point used for the analytic ex- 
trapolation: in our analysis we considered the amplitude given by the ex- 
pression (H), which was obtained by making the LSZ reduction of the kaon in 
the S-matrix element (|^). Alternatively, by reducing first one pion instead 
of the kaon, we obtain the expression 

Ai = ^^{<k,)\r,Mmp)). (54) 
v2A;io 

By further reducing the i^'-meson, one obtains, instead of (^), the expression 



Ai{s,t) = -^=_ j dxe-^P^^(xo) (7r(A;2)|[r/.(0),r/^(x)]|0) , (55) 

which allows the analytic continuation with respect to the variables s = 
and t = {p — k2Y- It is easy to see that the spectral function can be written as 
in Eq. (|^, with the corresponding terms similar to (§|)-(0), excepted that the 
momentum conservation reads now = p in as and Pn = ^2 — p va. at- This 
implies that the significance of the various terms in the spectral functions are 
different in the two approaches: in Eq. (|^) the cm. energy which generates 
the intermediate states is yielded by the total momentum carried by the kaon 
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and the interaction hamiltonian, while in (^) this energy is provided by the 
unphysical mass of the pion. In the alternative approach mentioned above, 
the energy in the term as is yielded by the unphysical mass of the kaon, 
while in at it is provided by the momenta of the pion and of the hamiltonian. 
Despite this different interpretation of the matrix elements, it is easy to see 
that the Omnes representation is formally similar in the two approaches. 
The differences concern only the inelastic contributions, parametrized by the 
expansions in powers of the conformal mapping variables. 



4 Discussion 

The Omnes representations derived above generalize previous results ob- 
tained in the literature. In Refs. Q the authors write down an Omnes 
representation for the decay amplitude as a function of s at fixed t = m^, 
using the formal analogy with the scalar pion form factor. The decay ampli- 
tude is identified, up to a constant, to the Omnes factor fli{s) in Eq. (p7D, 
i.e. the inelastic singularities are neglected. The role of the inelastic FSI 



and ISI contributions in the dispersion relation was discussed in Ref. |T2 
assuming that the decay amplitude defined in Eq. (|l|) satisfies a Mandelstam 
representation. 



A similar line is followed in Ref. [|T^], where the weak Hamiltonian Ti 



in Eq. (|ip is identified with a field operator with the quantum numbers of 
the kaon, and is assumed to carry a nonzero momentum. In this approach 
the K ^ niT decay amplitude is obtained from a dispersion relation for the 
elastic processes KK tttt (s-channel) and nK ttK (t-channel). As- 
suming that only S and P waves contribute to both channels, the amplitude 
is written in [|I^ as a sum of functions depending on a single variable, s 
(t), respectively. The interplay between the weak and the strong dynamics 
is however not apparent in this treatment: the weak K ^ mi decay am- 
plitude is given formally by the amplitude of the strong scattering process 
KK — > TTTT, evaluated at an unphysical point. 

As we mentioned already, the dispersion relations for the K —>■ tttt decay 
require a continuation in the external momenta. Therefore, the interpretation 
of the dispersive variables was not very clear in the previous works. In ref. 
PI the variable s in the dispersion relation for the amplitude Aj{s, t) at fixed 
t = was identified with the momentum squared of an off-shell kaon. This 



raised the subsequent criticism [13|, which emphasized the ambiguity of the 
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ChPT calculations for off-shell operators. We mention also that in Ref. [TT 
the same variable s was identified with the mass squared of the kaon. 

In the present treatment, the LSZ reduction formula allows the analytic 
continuation of the amplitudes in the complex planes of the external mo- 
menta. No off-shell extrapolation of the operators is necessary and the mean- 
ing of the variables is clear, as discussed at the end of the previous section. If 
we use as starting point of the analytic continuation the matrix element (P, 
the variable s is defined as s = {ki + /c2)^, and in unphysical configurations it 
may be different from the kaon momentum. Moreover, in this case the vari- 
able t is equal to kf, and represents the mass squared of the external pion. 
On the other hand, if we use as starting point of the analytic continuation 
the expression we have s = p"^, i.e. it represents the mass squared of 
the external kaon, while t = {p — /c2)^ includes, in unphysical configurations, 
the momentum carried by the interaction hamiltonians. This interpretation 
allows us to incorporate in the dispersion relations, at least to a certain ex- 
tent, the results of the lattice simulations. We consider for illustration the 
first definition of dispersive variables. 

Most lattice calculations simulate matrix elements of the form (nlOjlK), 
related to the matrix elements of interest {nirlOjlK) by lowest order ChPT 
in the soft pion limit 0. In the limit fci — we have t = kf = and 
s = {ki + ^2)^ = rri^. Therefore, the lattice simulations of this type provide 
the value of the amplitude Aj{m'^, 0). 

Direct simulations of the matrix elements {7^7r\Oj\K) are done only for 
special configurations, for instance when the kaon and one pion are at rest. 
Assuming that the reduced pion is at rest, ki = 0, we have s = t + + 
2\/tET^, where E^^ is the energy of the moving pion. At present the lattice 
simulations are done using values of the mass larger than the physical 
one. To match this situation we take i = ■fh'^, which means that s = ih\ -\- 

-|- 2rfi7r-E7r- We assume therefore that Ai{s,t) is known approximately 
from the lattice calculation^. Actually, Aj{s,i) does not correspond exactly 
to the configuration which is simulated on the lattice, since in the dispersion 
relation the unreduced pion has the physical mass, k\ = m^. The analytic 
continuation with respect to k^ requires the reduction of the second pion in 
the relation (^, and is more complicated. In the present formalism, we can 
assume that the lattice situation is approached by using a dispersion relation 



■^In the alternative interpretation of the dispersive variables, discussed above, the cor- 
responding points are s = mj^ and t = + — 2mA' Et^- 
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symmetrized with respect to the interchange of t and u, as in Eq. (p^, where 
one pion has the physical mass and the other a different mass, rh.,^. This is 
still not identical with the lattice case, but we believe that even approximate 
numbers are welcome, since the dispersive approach includes correctly the 
FSI. Using a sufficient number of points (s, i), it is possible to determine the 
coefficients c^^-* and c^^ of the Taylor expansions in powers of the conformal 
variables z and w, and to calculate then the physical amplitude using the 
expression (^8]) for s = m\ and t = m^. 



5 Conclusions 

In the present paper we derived Omnes representations for the K tttt am- 
phtudes, which include elastic and inelastic contributions in both the direct 
and the crossed channels. We showed that the amplitude is decomposed as 
a sum of two functions, one depending only on s and the other depending 
only on t. This decomposition follows naturally from the LSZ formalism and 
hadronic unitarity and does not require additional assumptions. The elastic 
contributions are parametrized by Omnes factors according to Watson the- 
orem, and the inelastic singularities are accounted for by the technique of 
conformal mappings. The treatment based on LSZ formalism allows a clear 
interpretation of the dispersion relations and the meaning of the dispersive 
variables. The unknown coefficients c^^-* and c^f^ entering the parametriza- 
tion (H) of the amplitude can be determined, at least approximately, using 
information provided by lattice calculations at unphysical momenta. The 
numerical implementation of this program is a subject of a future work. We 
mention finally that the effects of isospin violation, which were discussed 
recently in Ref. can be incorporated in the dispersive treatment by a 

suitable modification of the unitarity relation. 
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